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Abstract 

We consider a classical problem of linear stability of convective rolls in 
a plane layer with stress-free horizontal boundaries near the onset of con- 
vection. The problem has been studied by a number of authors, who have 
shown that rolls of wave number k are unstable with respect to perturbations 
of different types, if some inequalities relating k and the Rayleigh number R 
are satisfied. The perturbations involve a large-scale mode. Certain asymp- 
totic dependencies between wave numbers of the mode and overcriticality are 
always assumed in the available proofs of instability. We analyse the stabil- 
ity analytically following the approach of Podvigina (2008) without making a 
priori assumptions concerning asymptotic relations between small parameters 
characterising the problem. Instability of rolls to short-scale modes is also 
considered. Therefore, our analytical results on stability to space-periodic 
perturbations are exhaustive; they allow to identify the areas in the {k, R) 
plane, where convective rolls are stable near the onset. The analytical results 
are compared with numerical solutions to the eigenvalue problem determining 
stability of rolls. 
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1 Introduction 



We consider Boussinesq convection in a plane horizontal layer heated from below 
with stress-free horizontal boundaries. For small Rayleigh numbers i?, i.e. for small 
temperature differences between the lower and upper boundaries, the fluid is not 
moving and the heat is transported by thermal diffusion only. As R exceeds the 
critical value i?*, the fluid motion sets in. The motion takes the form of rolls. We 
denote by kg the horizontal wave number of the mode becoming unstable the first. 
(By a horizontal wave number we understand the length of the horizontal component 
of the wave vector.) In this paper wc study analytically and numerically stability of 
rolls of a horizontal wave number close to kg for the Rayleigh number slightly above 

Instability of rolls in a convective layer was studied analytically by a number of 
authors. An instability of rolls specific for stress-free boundaries is known, which 
does not occur in a layer whose one or both horizontal boundaries arc rigid. Its 
presence relies on existence of a slowly decaying large-scale mode. Zippelius and 
Siggia (1982, 1983) were the first to study the instability of this kind. In the leading 
order the unstable mode is a sum of a large-scale mode and of two short-scale modes 
with wave vectors close to the one of the perturbed rolls. In the study of stability 
of rolls in a rotating convective layer it was called small angle instability (Cox and 
Matthews 2000). 

By deriving amplitude equations, Zippelius and Siggia obtained sufficient condi- 
tions for instability of rolls, in particular, they found that for P < 0.782 (P denoting 
the Prandtl number) no stable rolls existed near the onset. Their results were ques- 
tioned by Bussc and Bolton (1984), who found boundaries for instabilities of rolls 
by direct calculations of the unstable mode, and claimed that no stable rolls were 
present near the onset only for P < 0.543. Their results were confirmed by Bernoff 
(1994), who studied instability of rolls applying Ginzburg-Landau equations. 

The conflict between the results of Zippelius and Siggia (1982, 1983) and Busse 
and Bolton (1984) was resolved by Mielke (1997), who studied stability of rolls by 
means of Lyapunov-Schmidt reduction and showed that instability boundaries had 
been found in these papers for different unstable modes. The problem of stability 
of rolls involves five small parameters: two wave numbers of the large scale mode, 
overcriticality, the difference between the wave number of perturbed rolls and the 
critical wave number, and the growth rate (depending on the first four). Zippelius 
and Siggia postulated asymptotic relations between the parameters, different from 
those postulated by Busse and Bolton, and hence different instability modes were 
examined. 

However, in all these studies some asymptotic relations between the small pa- 
rameters of the problem were assumed, and thus stability only to selected types of 
perturbations was studied. Hence, only instability of rolls was proven (as it is dis- 
cussed in introduction and conclusion in Mielke (1997) and also section 3 of Bernoff 
(1994). The question, whether other unstable modes corresponding to other asymp- 
totic scalings exist, remains open. Another question asked by Busse and Bolton 
(1984) and Bernoff (1994) and not answered by previous studies is whether enough 
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terms of asymptotic expansions were taken into account. For four independent small 
parameters this is a hard question! 

In the present paper both difficulties are resolved. We do not assume any asymp- 
totic relations between the small parameters, hence instability with respect to all 
perturbations of the small-angle type is examined. An unstable mode is represented 
as a series in small parameters, with estimates for the remainder. (Estimates for 
omitted terms were not given before.) A condition for instability has the form 
of inequalities. We demonstrate that in these inequalities the omitted terms are 
asymptotically smaller than the terms retained in the analysis. 

We also study instability to perturbations of a different kind, which are in the 
leading order convective rolls with the horizontal wave number close to the criti- 
cal one, and a finite angle between wave vectors of the perturbed rolls and of the 
perturbation, which we call a finite angle instability. 

Furthermore, we show that stability to all considered perturbations implies sta- 
bility to a much wider class of doubly periodic in horizontal directions perturbations. 
We show, that the domain of the linearisation operator, acting on vector fields sat- 
isfying the assumed boundary conditions on the horizontal boundaries and doubly 
periodic in horizontal directions, splits into a direct sum of invariant subspaces. 
Thus investigation of instability is reduced to detection of instability modes in these 
invariant subspaces. Any instability mode in such a subspace is either of a small 
angle or finite angle type. Therefore, our study of instabilities of rolls with respect 
to perturbations which are doubly periodic in some directions in the {x, y) plane is 
complete: we demonstrate stability of rolls for which instability is not detected in 
our analysis. This is also a novel feature of our investigation: all previous papers 
focused exclusively on instability. 

We estimate analytically the asymptotics of the most unstable mode and its 
growth rate on different parts of the {k, R) plane and calculate the instability bound- 
aries estimating the orders of the neglected terms (which was not yet done in liter- 
ature) . 

We are using notation and some results of Podvigina (2008), referred to as 
OP2008, where instability of fiows in a rotating convective layer with stress- free 
boundaries was studied. 

Stability of rolls in a layer with stress-free boundaries was studied by direct 
computations of dominant eigcnvahics by Bolton and Bussc (1985). Their numerical 
results agree well with the theoretical predictions of Busse and Bolton (1984), in 
particular, they found that for P = 0.71 rolls are stable near the onset. However, 
Mielke (1997) claimed that stable rolls near the onset for P < 0.782 are absent; 
he did not comment on the disagreement with the numerical results of Bolton and 
Busse (1985). In contrast to (Bolton and Busse 1985), our computations of stability 
modes indicate that rolls are unstable at the onset for P < 0.782, albeit in a small 
neighbourhood of the point {kg, R^) on the {k, R) plane. For P decreasing below 
0.782, the right boundary of the area of stable rolls slowly moves to the left, away 
from the point {kg, R^). 
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2 The onset of convection 



V X (V X v) + PV V + PROe, - Vp, 



(2) 



(3) 



Boussinesq thermal convection satisfies the Navier-Stokes equation 
dt 

the incompressibility condition 

V- v = 

and the heat transfer equation 

f)f) 

^ = -(v-V)^ + t;. + V2^ 
at 

where v is the flow velocity and 6 is the deviation of the flow temperature from 
the linear profile. R and P are dimensionless parameters, the Rayleigh and Prandtl 
numbers, respectively. Stress-free horizontal boundaries held at fixed temperature 
are assumed: 

dvx dvii ^ , , 

^ = ^ = t^. = 0, ^ = atz = 0,l. 4 
oz oz 

The trivial solution (v, 9) = (0, 0) describing pure thermal conduction loses 
stability to perturbations of wave number k at R = Rc{k), where 

R^{k) = a^k-^, a = k'^ + vr^ (5) 

The critical horizontal wave number ks for the onset of convection is tt /V2, the 
respective critical value R^ = 277r^/4. 
We employ four- dimensional vectors 



W = (W*^°", W^^'^P) = (v,e). 
For a Rayleigh number slightly above the critical value, 

R = Rq + 

a solution to (II])- (Hj) can be represented as a power series 

oo 

The first two terms of the solution representing rolls are 

/ —Trk^^cosTTZsinkx \ 


sin TTZ cos kx 
\ sin irz cos kx J 

[Ui is an eigenvector of the linearisation of dl])-® ), 

/ \ 



(6) 
(7) 



Ui 



(9) 








\ -(Svra)-^ sin27r5; / 



where 



(10) 
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3 Invariant subspace 



We study stability of rolls of wave number k, which is close to the critical one; 

a = k-ks (12) 

is thus a small parameter. 

Stability of (IHl) is controlled by eigenvalues of the linear operator L, a linearisation 
of ([I])-© near the steady state. The operator can be expanded in a series 



j=0 

Here the first term is 

Lo(v, 6) = {PV\ + PRcOe, - Vp, v, + 
We consider the eigenvalue problem 

LW = AW. 



(13) 



(14) 



(15) 



In OP2008 a three-dimensional invariant subspace of L was considered, and the 
problem of stability of rolls was reduced to the analysis of eigenvalues of L in this 
subspace. Denote by W^, j = 1, 2, 3, a basis in this subspace and by A the matrix 
Aij of the restriction of L on the subspace: 



LWj = J2AijWi, J = 1,2,3. 

i=l 

We expand the basis and the matrix in power series in e: 

oo 



(16) 



(17) 



1=0 

oo 



A: 



1=0 



The coefficients in the series (fTTl) and (ITSl) depend on 6^ and 6y, which are small 
perturbations of the horizontal component of the wave vector (/c,0,7r): 6^ k and 
6y <^ k. (The two small parameters are involved in the definition of the three- 
dimensional invariant subspace of L, see OP2008.) 
Vector fields W^ q are eigenfunctions of Lq: 



(19) 



the leading order approximations of Wj^o in and 6y were derived in OP2008: 

/ -6y sm{6xX + 5yy) \ 
5x sin(5^x + 6yy) 


V / 



Wi,o 



(20) 
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/ 



ns^kj^^ COS nz sm{{k + 6x)x + 6yy) \ 
-Tx5yk^^ cosTiz sin((A; + 5x)x + 5yy) 



W2,o 



sin TTz cos((/c + 5x)x + 5yy) 
simx z cos{{k + 5 x)x + 5 yy) J 



+0{{kl-k^)\a{kl-e)), (21) 



v 



--Ks_k_^ cosT[zsm.[{k — 5x)x — 5yy) \ 
n5ykZ^ cos Txz sin((A; — 5x)x — Syy) 



W 



3,0 



sin nz cos((A; — 6x)x — 6yy) 
aZ^ siiiTcz cos{{k — Sx)x — 5yy) J 



+ Oi{k^-ky,a{e_-k^ 



)) (22) 



v 



(the expression for Wi^o is exact). Here it is denoted k± = ((/c ± 6xY + SyY^^, 
s± = {k ± Sx)k^^ and a± = k"^ + n^. Wi q is a large-scale horizontal mode. 

The paper is mainly concerned with the eigenvalue problem ( ITSjl in the three- 
dimensional invariant subspace, discussed above, where the eigenfunction W and 
the operator L take the forms (ITTll and ( ITSll . respectively. The case of an unstable 
mode from this subspace is called small angle instability, because the angles between 
the wave vector {k, 0, vr) of perturbed rolls and the wave vectors {k + 6x, Sy, vr) of 
short-scale components of perturbation (W2,o and W3 q in the leading order) are 
small. In section 5 we also study stability with respect to perturbations of the form 
of rolls at an angle ^ to the perturbed ones, where ^ is finite (and not asymptotically 
small). In the remaining part of the section it is shown that any growing mode of 
L belongs to one of these two classes. 

Consider the space J^{6x,Sy) of 4-component vector fields of the form ([6]), such 
that Wfl°™ and W*^"^p 

(z) are linear combinations of harmonics with wave vectors {rriik + 5x,Sy,'m27T) or 
{mik — 6x, —6y,m2n), where mi and m2 are integer, 
{a) are symmetric about the vertical axis, i.e. 

(v(x, y, z), e{x, y, z)) = {-Vxi-x, -y, z), -Vy{-x, -y, z),v^{-x, -y, z), e{-x, -y, z)), 

{Hi) satisfy the boundary conditions (jl]), 
{iv) W^°™ is divergence-free. 

It was shown in OP2008 that J-" is L-invariant for any (not necessarily small) Sx and 

6y. 

Let the subspace ^{Sx, Sy) be defined like we have defined J-'{Sx, Sy), but omitting 
the condition {ii). Such subspaces also are L-invariant. We consider perturbations 
which are doubly periodic on the {x, y) plane. Let the domain of L be comprised 
of 4-component vector fields, for which [Hi) and (iv) are satisfied, and which have 
the same double periodicity on the {x, y) plane. The domain can be split into a 
sum of invariant spaces J^{6x,Sy); hence we can assume that a mode belongs to 
such an invariant subspace. Any mode can be represented as a sum of a symmet- 
ric and an antisymmetric vector field; each of these vector fields is itself a mode, 
since the subspaces of symmetric and an antisymmetric vector fields are L-invariant. 
Moreover, in a coordinate system with the origin shifted by half a period, ir/Sy, in 
the direction (this shift does not affect rolls, since they are independent of the 
y coordinate) the antisymmetric modes become symmetric. In the case Sy = 0, if 
integer m and n such that 6x/k = (2m + l)/(2r2) exist, the shift of the origin by 
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/ = 27^71 /k = 7r(3m + l)/6x along the x— direction turns a symmetric mode into an 
antisymmetric one. For a given k and Sx the ratio (2m + l)/(2n) can be arbitrary 
close to 6x/k. Consequently, without any loss of generality we consider henceforth 
only modes belonging to J-". 

Eigenvalues of f ll3p are perturbations of the ones of Lq. Positive or slightly 
negative eigenvalues of Lq are associated with eigenvectors, whose wave vectors are 
either {Sx,Sy,0) (large-scale mode) or (/ccos^, /csin^, vr) with k close to kg (see, e.g., 
discussion in Bernoff 1994). The latter eigenvectors are (Q rotated by the angle ^ 
about a vertical axis; we denote them by U(/c,^). Consider an eigenvector W of 
L which belongs to some J-". The vector field Wq = limg^oW can be one of the 
following: 

(a) Wq = ai Wi + a2W2,o + c^aWa q (with at least one aj ^ and Wj^o defined by 

dloD-dliD ); ' ' ' 

(b) Wo = U(A;, with k close to fc„ ^ = 0(1) and i ± 27r/3 = 0(1); 

(c) Wo = aiU(A;,^i) + a2U(A;,^2) with A; close to A;,, 6-27r/3 = o(l) and 6 + 27r/3 = 
o(l). 

In case (a) small-angle instability takes place studied in section 4, in cases (b) and 
(c) finite angle instability considered in section 5. Therefore, we examine all types 
of possibly growing perturbations. 

4 Stability of rolls: analytical results 

It is shown in Appendix B that in order to study stability of rolls with respect to 
perturbations from the subspace constructed above, it suffices to check, whether (for 
a given P, e and a) there exist such and by that the following inequalities are 
satisfied: 

det ^ > (23) 

or 

5(^)tr^-det^>0 (24) 

where 

S{A) = A11A22 - A12A21 + - + A22A33 - A23A32 

is the sum of the three second order minors. The matrix A has an eigenvalue with 
a positive real part, if and only if at least one of the inequalities is satisfied for some 
(5's. 

It is calculated in Appendix A, that in the leading order 

det A = Dq + e^Di + e^D2 + aDs + a^D^ + ae^D^, (25) 
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where 



Do 



D2 



D 



'3 



{51 + 5^y)-\d,r5l + d,25t5l + d^^Sld'^ + ^14^,^), 



S{A)iiA- det ^ = + £^^1 + e^E2 + a^Eg + ^£^^4, 



(26) 



where 



^3 = es^^ 

i?4 = 645^52(^2 + ^2)-l_ 



The values of dij and Cjj are given by expressions (1671) and (l69l) in Appendix A. 

To investigate stabihty of rolls we consider exhaustively different asymptotic 
relations between a and e and different signs of a. We also consider two limit values 
of the Prandtl number. 

Examples of areas on the (fc, R) plane where rolls are stable are shown on Fig. 1 
for several values of P. The area of stable rolls found numerically (see section 8) 
is shaded; the instability boundaries determined analytically are shown by lines. 
We use the standard notation (Busse and Bolton 1984; Bernoff 1994; Mielke 1997; 
Getling 1998) for the instability modes and respective instability boundaries. For 
a skew-varicose (SV) mode the associated eigenvalue is real, and for an oscillatory 
skew-varicose (OSV) mode the associated eigenvalues are complex; both modes exist 
for 6x ~ 5y. A zigzag (ZZ) mode emerges for 5x = and the associated eigenvalue is 
real. The mode responsible for the instability for P < 0.782 exists for 6y ^ 6^ ^ 
and the associated eigenvalue is real; it is also called a skew- varicose mode (see, e.g., 
Mielke 1997). To distinguish this skew- varicose mode from the SV mode, we label 
the former SV2. 



4.1 The case <C 

Let 6l ^ 5y and 5l <^ e^. Then in the leading order 

det^ = £^^21 5,^. 



(27) 
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d2i is positive for all P (see (!67|) ). hence det^ > as well, implying that for 
the assumed (5's the matrix has a positive real eigenvalue and the rolls suffer from 
monotonous instability (it can be shown that oscillatory instability does not emerge 
in this case). 

4.2 The case a> 

Assume again 5^ ^ 5^ and 51 <^ e'^. In the leading order the determinant now 
includes a— dependent terms: 

det A = e'^d2i6l + aD^ + a^D^ + ae^D^. 

However, for a > the terms involving a are positive. The first term is positive, as 
discussed in the previous subsection, hence rolls are also monotonously unstable in 
this case. 

4.3 The case a < 0, P < Pi « 0.782, ~ 

Denote by Pi the Prandtl number which is the solution to the equation 

i.e. 

-p2 + 2P + 2 = 2(2P2(P+ 1))^/^ (28) 

The solution is 

Pi = (3 - 2V2)(1 + ^7 + 4^2) ^ 0.782 

(cf. Zippehus and Siggia 1982). If 5l < 5^, 5l > 5y and ^ e^, the sum of 
asymptotically largest terms in (1251) is 

d,25l5l + e^di:,5l5l + e'd225l (29) 

which is positive for P < Pi. Hence, for P < Pi and the assumed a's rolls are 
unstable. 

4.4 The case a < 0, P > Pi, ~ 

Re- write (12^ as det^ = Vi +'D2, where 

Vi=Do + e^Di + e%2Sy{Sl + 6l)-' + aD^ + a^D^ + ae''d^25l (30) 

V2 = e\eH2i + c^d,i)5l5l{5l + 5j)^\ (31) 
As proved in Appendix C, Pi < for P > Pi and the assumed a's. P2 > if 

SV: e^>-a/i, /i = ^ = ^ A (32) 



9 



For ^ d^.^ 6y V2 is asymptotically larger than Vi, hence (132|) yields a boundary 
for monotonous instability of rolls. 

Represent (l26l) as S{A)tYA — detA = 81 + 82, where 

8i = Eo + e^Ei + a^Es (33) 

82 = e\e^E2 + aE4). (34) 
£^1 < for P > Pi and assumed a's. £^2 > if 

^ / \ 1 /2 

622 + > 2(621623) ' , 

108(P + l)27r2A; 



(P + 3)(3P2 + 2P + 2) + 3P2(P + l)i/2(p + 5)1/2 ■ 

For ^ 5l and ^ 5^, is asymptotically larger than 81, hence (|35|) is a 
condition for instability. 

4.5 The case a < 0, ;» 

Assuming, as above, e:^ ^ 52 g^j^^^ ^2 ^ ^2^ g^^^^ ^j-^^^^ e^Q,^^ > is asymptotically 
the largest term in (l34l) . Hence, rolls are unstable. 



4.6 The case a < 0, ~ and large P 



For P — 7- 00, /i in fl32l) has a finite limit, and /2 in fl35|) vanishes. Hence, for 
P sufficiently large the instability under the condition ( 135|) can compete with the 
instability occurring for a ~ e. As shown in Appendix D, for large P, det ^ > 0, if 

q7r2p2 

ZZ: ^^</3a^ /s = (36) 

4.7 The case a < 0, ~ and P sUghtly smaller than Pi 

Mielke (1997) showed that for P < Pi rolls near the onset are always unstable. More 
precisely, the following has been proved: for such P there exists a neighbourhood 
of the point [kg, R^) in the {k, R) plane, such that for a given R rolls of horizontal 
wave number k (where they exist) are unstable. The question, how the area where 
rolls are stable is modified, as P becomes smaller than Pi, has not been addressed in 
literature. We show below that the area of stable rolls does not disappear abruptly 
and it still exists near the onset, but its boundary does not include the point {ks, P*) 
(cf. figs. la,b and figs. Ic-f). As P decreases, the area of stable rolls moves away 
from this point, because the SV2 boundary moves to the left. 

In the search of the horizontal scale ratios for which (!25ll is positive for P < Pi, 
it was assumed in section 4.3 that 5^ ^ 5^ ^ 5^ and then in the leading order det A 
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is given by (l29l) . For P slightly smaller than Pi the expression (!29l) can be of the 
same order as other terms in (!25|) not far from the onset. Consider a new small 
parameter P = Pi — P. The maximum of (!29|) is admitted for 

« = - ie\ (37) 



and the maximum is equal to 

8l7r^P(P + ip " 

where 



_ 4(-2(-2Pi + 2)(-Pf + 2Pi + 2) + 24Pf + 16Pi) ^ ^ . 

8l7r4A(A + l)2 ^-^'^^ • 

For 5^. defined by (1371) in the leading order the determinant is 

e^{e^{e'^{dQi + do2)g^ + £^^12?^ + £^^21? + a^sig) + ^^45^/3 + (c^oag + c^m)^^) (38) 

which is a cubic polynomial of 5y. Its maximum is admitted, when 

54 ^ ^^^6^ (39) 
y 3{do3q + du) ' 

and the maximum of (1381) is 
where 

S3 = d^iq ^ 0.0299, si = s^\{doi + c/o2)g^ + d^q^ + c/aig) ~ 0.00279, 
2,1/ -4 \ 

The maximum is positive for 

a > -e'^si - efl'^S2. (40) 



However, the boundary of the SV2 instability defined by f l40|) turns out to be in 
a poor agreement with the numerical results discussed in section 8. For P = 0.6 
and P = 0.7 (figs. la,b) the right boundary ( HOl) of the area of stable rolls is shifted 
far to the left compared to the computed one. The SV2 boundary defined by the 
condition max^^^^ det^ = with all leading terms in (!66l) retained is still shifted 
too far to the left. The asymptotics fails because the values (1371) and (l39l) of bx and 

respectively, are of the order of 0.1 for the considered overcriticalities ~ 1, 
while the asymptotic analysis is applicable for infinitesimally small bx and by. (For 
example, as we have found numerically, for P = 0.7 the intersection of the SV and 
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SV2 boundaries is at a = —0.00172 and = 0.583. The respective values (!37|) and 
(139!) are 6x = 0.126 and Sy = 0.413. For the SV2 dominant mode on the stabihty 
boundary the computed values are also large, 6^ = 0.124 and 6y = 0.332, which 
surprisingly do not differ much from the values obtained analytically.) 

Consequently, we follow an alternative approach and assume that the SV2 in- 
stability boundary can be described as an equation, where it suffices to retain two 
first non- vanishing terms in the Taylor expansion in e and (3. The condition for the 
instability thus takes the form 

SV2 : a> e^hi + £/3/i2, (41) 

for some coefficients hi and /i2, which can be determined numerically. The intersec- 
tions of the SV2 boundary with the SV and OSV boundaries have been computed 
for P = 0.6 and P = 0.7 (see figs. la,b). The minimum (over the coefficients hi) of 
the maximum (over the four points of intersection) relative error is equal to 0.14, it 
is admitted for 

hi = -0.0012 and /la = -0.018 . (42) 

By the relative error we understand the ratio | (etc— a*) /«c| , where ac is the computed 
value (see section 8) and at is found from (HTll at the points of intersection. We have 
also computed several points on the SV2 boundary in the regions of other instabilities 
of rolls, employing the fact that the respective (local) maximum of A over 6x and 6y 
is admitted for 6y much larger than for the other instabilities. The computed values 
agree well with (HTl) . P2l) (see figs. la,b). Fitting of si and S2 in PUI) yields a much 
higher (about 0.5) minimum of the maximum over the four points relative error. 



5 The finite angle case 

To analyse stability of rolls of horizontal wave number k to rolls of wave number 
kp, which are rotated by angle ^ with finite ^, we use center manifold reduction 
with the center eigenspace spanned by eigenvectors ([9]) with wave vectors {k, 0, vr) 
and (/Cp cos.^, fcp sin^, vr). We perform the reduction like in (Podvigina and Ashwin 
2007). Here only results of calculations are presented. Periodicity in horizontal 
directions of the considered rolls implies that periodicity cells are parallelograms 
(and not squares as ibid). 

Restricted to the two-dimensional (C^) center manifold, the system has the form 

(43) 

Z2 = X2Z2 + 2;2(^3klP + ^4^2^), 

where zi and Z2 are coordinates in the center manifold, along the directions (A;, 0, vr) 
and (fcp cos^, fcp sin^, tt), respectively. The reduction is performed for R = Rc{k). 
We are interested in kp close to k (otherwise A2 is of the order of one and negative 
and thus the rolls (zi, 0) are stable near the onset). 
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For e defined by ([7]) and kp close to the coefficients of linear terms in (143|) are 

Ai = P{P +l)-^k^a-^e\ 

A2 = P(P + (^A;2a-2e2 _ 4(A;p _ kf - 8a{kp - A;)) . 

For small kp — k the differences Ai — A4 and A2 — A3 are small, coefficients of cubic 
terms in f H3|) are 

0.125P 

~ = -(pTT) • 

0.125P a(l - cos^^) .,2ag+ + 2Pa2 2aq^ + 2Pa? 



A3 ~ = r —, (1-COsO . + (l+COSn 

^ P + 1 3P+1 ^' PA, ^ ^' 



(P + 1) 3(P+1) V' PA+ ' PA_ 



■ (^(1 - cosO + (1 + cosO 



2(P+1)V^ A+ ■ V ■ ^/ ^_ 

where 

/i± = 2A;^(1 ± cosO, q± = 47r^ + h±, A± = - Rh±. 
The amplitude of emerging rolls is 

and four eigenvalues of fll5|) linearised around the steady state are 

-2Ai, \2 + As\zi\^, 0, 0, (44) 
hence the instability condition is 

e^ea-\l - A^V^i) - 4{k - kpf - 8a{k - kp) > 0, 
i.e. instability occurs if 



AAia'^ 



^3 ) 

where by A™'*^ we have denoted the maximum of ^3 in ^. 

For a finite P the instability boundary is below the boundary defined by (l35l) . 
For large P, the limits of Ai and ^3 are finite, hence /s < /3, and finite-angle 
instabilities do not affect the area of stability of rolls. The instability with respect 
to rolls rotated by ^ = 7r/2 is called the cross- roll instability. Note that the maximum 
of A3 can be admitted for a ^ 7^ 7r/2, but we do not consider here the problem of 
maximisation of ^3 in ^. 

For ^ close to 27r/3 the center eigenspace also involves rolls with the direction 
of the axes rotated by — 27r/3. The system restricted to the three-dimensional (C^) 
center manifold is 

ii = Ai2;i + + A2\z2\'^ + Bi\z3\'^), 

Z2 = X2Z2 + 22(^3klP + ^4^2^ + B2\Z3\'^), (45) 
Z3 = X2Z3 + Z3{A5\zi\^ + A6\Z2\^ + B3\z3\^). 

However, the eigenvalues determining stability of rolls are (jSl), examined above. 
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6 Growth rates 



In this section we find orders of growtli rates of the dominant unstable modes. If 
entries of the matrix A are of different orders, it is possible to calculate dominant 
eigenvalues, like it was in the case of rotating layer in OP2008. In the present 
problem often almost all entries of A turn out to have the same asymptotics, and 
only orders of growth rates can be determined. Also we find orders of coefficients 
ttj, j =1,2,3, of the most unstable mode 

W = aiWi + asWa + 03 Wg. (46) 

Unstable modes (or instabilities) can be roughly categorised into five different 
type^: 

SV : (5a; ~ Sy 

SV2 : 6x < 6y 

OSV : 4 ~ Sy 

ZZ : 6x = 

E - 1 : 2k6x + 51 = -2ka 

Similarly to section 4, we consider different asymptotic relations between a and 
e. Our findings are summarised in Table 1, where 

_ (p2_2P-2)2-8P2(p+i) _ 4(2A;)i/2 . _ 4P t _ ^ 
~ 97r2(P + l)2p ' ^2 - 3vr(P+ 1)1/2' « " (p + i)' " ^^^^ 

^E-1 stands for Eckhaus-like instability. Maximisation of the growth rate in Sx and 6y yields 
only the horizontal wave number of the most unstable mode, see OP2008 and section 6.5. The 
conventional Eckhaus instability is a particular case of the E-1 instability for Sy = 0. 
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Table 1. Possibly dominant instability modes for various asymptotic relations be- 
tween a and e and values of P. The last column presents eigenvalues, when they 
can be calculated, or their orders of magnitude otherwise. (Hence, often it remains 
unclear which mode is dominant.) 



Relations 


Conditions for 


Type of 


6x and 6y 


Eigenvalues 


between a and e 


existence 


the mode 








none 


SV 


5x ^ 5y e 






P < Pi 


SV2 


5x ^ ^ £ 


A = 




> -fia 


SV 


^ Sy e 


A~£2 




P < Pi 


SV2 


5x < 5y, 5x ~ e 


A = 




< -/2a 


osv 


Sx ^ 5y e 


Re(A) ~ £2 


< a2 < £4/3 


a > 


SV 




A = 


or 


a < 0, 


osv 


51 ~ 51 ~ ea^l^ 


Re(A) ~ 


«2 ^ £4/3 


£^ < — /2« 








0? ~ 


a < 0, 

£2 > -/2a, 


zz 


(5a; = 0, 5y ~ e 


A = -^4^' + e3«' 


0? » 


none 


E-1 


2A;5^ + 51 = -2ka 


A = e3«' 



6.1 The case <C 

In this case there exists a growing mode 

SV : A ~ e^, for 5l ~ 5^ ~ 

since it can be easily shown that det ^ > for some 

5l^5lr^ e\ (47) 

If Wl\ holds, tr ^ ~ e^, SA ~ and det A ~ e^, implying that eigenvalues are ~ e^. 
For the assumed dependence of 5's on e, after the change of variables W2 — )■ eW2 
all coefficients (except for ^32) become of the same order in e, implying that for the 
associated eigenmode 01/02 ~ 03/02 ~ 

For P < Pi the expression ( l66l) can be positive also, if 5^ 5^ ~ e^. For 
5l + 5y^ e, the matrix (1651) has an eigenvalue close to — P(5^ + 5^). In the leading 
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order the associated eigenvector is Wi + ^2W2 + ^sWa, where ^2 = ^21/(^11 — ^22) 
and ^3 = A^i/^Aii — A33). Two remaining eigenvalues are eigenvalues of the matrix 



^22 " ^2^12 ^23 " ^2^13 
^32 ~ ^3^12 ^33 ~ ^3^13 



(4^ 



In the leading order they are 



y {51 + 5if 



± 

where 



(49) 

1/2 



Calculating their maxima in 5^ and 5y we find the most unstable mode: 



SV2 : \ = e 



2vP^-2P-2)2-8P2(p+i) 
97r2(P + l)2p ■ 



for 



,8(P + 1)2- (p2-2P-2)^ 



-, 5y > 5x 



18(P+1)P2 

These relations between 5's, e and A imply that the associated eigenvector (oi, 02, 03) 
of the matrix A f l65|) has components with the asymptotics ai/a2 ~ ^^"^j^^ ^'^d 
03/02 ~ e. 



6.2 The case ~ 

Dominant eigenvalues and asymptotic relations between the coefficients oi, 02, 03 for 
the eigenmodes SV and SV2 are the same as above. The SV mode is growing if f l32p 
holds, and SV2 if P < Pi. 

A growing oscillatory mode can exist, if holds true. As discussed in Ap- 
pendix B, condition f l2^ does not guarantee its existence. However, if such mode 
exists for all ^ e:^, a < 0, (see section 6.4), by continuity it exists for some 

~ e^. The maximal growth rate of the mode 

OSV : Re(A) ~ 

is admitted for 

X2 x2 ^2 

By the same arguments as for the SV mode, for the OSV eigenmode with the 
maximal growth rate and ai/a2 ~ 0.3/0-2 ~ e. 
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6.3 The case e^/^ > > a > q 

We employ the same change of variables as in section 6.1, and maximisation in 5^ 
and 5y yields the maximal growth rate for SV modes 

SV : A = = 3.r(P + l)i/2 ^"'^^ = « (50) 

The eigenmode coefficients have the asymptotics ai/a2 ~ ea'^^"^ and 03/02 ~ 
e'^{a^^'^5x)~^ ■ The SV2 mode has the growth rate 0(£^,a^), which is asymptoti- 
cally smaller than (!50|) . 

6.4 The case e^l^ > > e^, a < 

For the dominant oscillatory mode the maximal growth rate and (^x^^y) for which 
it is admitted are: 

OSV : Re(A) ~ at 5^ ~ 5j ~ ea^/^. (51) 

This can be obtained by the following arguments. Assume 

5l^5l<^ ea^'\ (52) 

Consider the cubic equation det(AI— ^) = 0. The assumption fl52l) implies, by virtue 
of the standard formulae for roots of cubic equations, existence of complex roots 
with a positive ReA ~ (det AY^^ ~ e'^a5l5y{5l + 5y)~^. Since this holds true for any 
{5x, 5y) satisfying fl52|) . this relation remains true for a 5^ ~ 5^ ~ ea^^'^. The last two 
asymptotic relations imply that the associated eigenvector (01,02,03) of the matrix 
A (jnS]) has components with the asymptotics 01/02 ~ a^/^ and 03/02 ~ a^/^e"^/^. 
For the SV2 mode the growth rate is ~ or ~ a^, i.e. it is asymptotically smaller 
than §B>. 

6.5 The case > ^4/3 

Maximisation of fHUl) in 6^ and 6y yields that the maximal growth rate 

E - 1 : A = AC^k^a^ (53) 

is admitted for 

(kl - A;2) = -2ak; 

the associated eigenvectors are either W2 or W3. (Note that growth rates are 
asymptotically smaller than fl53|) . unless + 6y ^ e; if this asymptotic relation 
is satisfied, employed in maximisation is valid, see section 6.1.) Alternatively, 
(I53p can be obtained directly from flMl) . For the OSV mode Re(A) ~ e^^^, which is 
asymptotically smaller, than (!53|) . 
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6.6 The case ~ e^/^ 

The maximal growth rate is O(a^) ~ 0(e^/^). For a > it is given by f l50|) . if 

or by otherwise. For a < it is either (IFIl) . or (155]) . 

6.7 The case ~ e^, a < and large P 

As noted in section 4.7, in the hmit P — t- oo the ZZ instabihty with 5x = ^ competes 
with the OSV instabihty and becomes of importance near the onset. If = and 
5y ^ £, the eigenvalues of the matrix A are 

where 

f2 , r4 



For large P, A3 > A2, maximisation of A3 in 5y yields the maximal growth rate 



Amax = -£^&^T7^ + 4C4/C^al 



4PA;2 

The associated eigenvector has asymptotics 01/03 ~ 1 and 02 = 0. 

7 Asymptotics of neglected terms in equations for 
stability boundaries 

Expressions (!32l) . (!35l) and (!36l) determining stability of rolls are only asymptotically 
correct. In this section we estimate the asymptotic order of errors in calculation of 
boundaries, relying on the known orders of the remainder terms in (1661) and (1681) . 

In the course of derivation of an equation defining the SV instability boundary, 
detv4 has been expressed in section 4.4 as a sum of Vi fl30|) and T'2 (I3T1) . where Vi 
is negative and involves terms 0(5^ + e^5'^ + e^5y5~'^) (here 5"^ = 5l + 5y) and V2 can 
be positive and involves terms 0{e^5l5y5'''^). The inequality ( l32l) is a restatement 
of the condition V2 > Near the boundary £^ 3> 5^ ^ 5y must be satisfied so 
that the sum (13T|) were positive. Under this condition Di is asymptotically smaller 
than T>2 and hence asymptotic corrections to T>i do not affect the boundary. Upon 
reintroduction of the terms omitted in (166!) . that are not asymptotically smaller than 
T^i^ dSID becomes 

V2 = e^{e\d2i + 0{e\ a)) + a{d^^ + 0{e\ a)))5l5l{5l + 
and thus the equation for the boundary takes the form 

SV: £^ > -a/i + 0(a^). (54) 
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Since at the boundary the factor in front of 5'^5y5 ^ vanishes, it can be shown using 
this analysis that the condition 

max det ^ = 0, 

defining the boundary, imphes 5^ = 5y = 0. 

The OSV instabihty boundary fl35p has been found from the condition that £2 
vanishes. By the same arguments as above, near the boundary ^ 5^ ~ 5y 
must be satisfied for the sum of £1 and £2 to be positive and hence asymptotic 
corrections to £1 again do not affect the boundary. With the omitted terms of fl68l) 
reintroduced, (IMD becomes 



£^=e^(^6l+6l)-\e\e2i+0{e\a))6l+{e^e22+ae,+0{e\a'))5l^^^^ 

This expression results in the equation for the instability boundary in the form 

OSV: <-af2 + 0{a^). (55) 

Again, it can be shown that at the boundary 5x = 5y = Q. 

In Appendix D the ZZ instability boundary is calculated from the condition that 

max det ^ = 0, 

5y 

where det^ is given by (1751) . Since on the boundary ~ e^, the equation with the 
omitted terms reintroduced takes the form 

max(-5^(2e2c3 + COie'h''^^ + PC,0 + 0(5J^ a^^J, ^6^2)^ ^ ^ (53) 
The maximum is attained for 

5l = 2ak + O(a^) 
and the condition for the instability is thus 

ZZ: < fsa^ + 0{a^). 



8 Stability of rolls: numerical results 

To examine stability of rolls of wave number k, we solve numerically (with an adapted 
version of the code of Zheligovsky 1993) the problem f[T^ for the eigenfunction 



f „..l pimkx 



w 



m=M n=N 

E E 

m=-M n=0 



COS Tinz 



\ 



w^„e""'^^ cos nnz 
w'^^e^"^^^ sin nnz 
w^^e^"^^^ sin nnz 



(57) 
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In computations, the cut off of tlie series at = M = 15 suffices (the spectrum of 
the solution in the Fourier space decays by at least 12 orders of magnitude). Location 
of the maximum of Re(A) in 5x and 5y has been determined with the precision of 
10~^ (or 2 ■ 10~^, if 6x and 5y are below 0.01) which allows us to find correctly at 
least two significant digits of A. 

The dominant eigenvalues of (IT^ . (|H7|1 and the values of 5^ and 5y where the 
maximum is admitted are shown on fig. 2 for P = 2 and k = 2.15 (thin vertical 
line on fig. Ic). In the interval 658.5 < R < 658.8, i.e. for a small overcriticality, 
the Eckhaus mode with = is dominant, the values of 6x and A are close to the 
ones given in the Table (according to the Table, A = 0.0136 and Sx = 0.071). In the 
interval 658.9 < R < 670 the dominant eigenvalues are complex, they are associated 
with the OSV eigenmode. The change of the type of the dominant mode implies 
a discontinuity of 5's. In the interval 658.9 < R < 661, Re(A) depends linearly 
on e = {R — RcY^"^, in agreement with the Table (in fact, for smaller R, where the 
instability is subdominant, this asymptotics for the eigenvalue of the OSV mode was 
also confirmed numerically), as predicted in section 6. For higher R the dependence 
is different, because, as noted in section 7, near the SV and OSV boundaries Sx 
and 6y become asymptotically smaller than e, while in section 6 we have assumed 
6x = 0(e) and 6y = 0(e). In agreement with section 7, we observe that Re(A), 6x and 
6y vanish at -Rosy and -Rsv, where -Rosy and i?sY denote the critical values of R for 
the OSV and SV instabilities. In the interval (-Rosy; -Rsy); where rolls are stable, the 
maximal growth rate is zero, admitted for Sx = Sy = 0. For R > i?sY the SV mode 
is dominant. Near -Rosy and -Rsy; Re(A) ~ (-Rosy ~ -R)^ and Re(A) ~ (-Rsy ~ -R)^; 
respectively (while no power law asymptotics has been found for 6x and 6y, except 
for Sy is almost linear near -Rsy)- Consequently, the SV and OSV boundaries are 
found by linear extrapolation of (Re(A))^/^ through two computed points close to 
the boundary. Near the ZZ and SV2 boundaries A depends on -R linearly, and we 
find the instability boundary by linear interpolation. 

The areas of stable rolls found numerically are shown on fig. 1 for several values of 
P. The difference between the SV, OSV and ZZ boundaries predicted theoretically 
and found numerically agrees with the estimations of the remainder terms obtained 
in section 7. For small a and e, the theoretical and numerical boundaries visually 
coincide, and the discrepancy remains small on increasing a. The area of stable 
rolls found numerically is shifted up compared to the one determined analytically, 
indicating that the contribution of the omitted terms is positive for OSV and ZZ 
instabilities, and negative for the SV instability. 

In view of the good agreement of the analytical and numerical results for these 
three boundaries, the disagreement for the SV2 boundary is surprising. A possible 
explanation is that the omitted in fl66l) terms involving Sx (which are of no impor- 
tance for the SV and ZZ instabilities for which Sx = - e.g., Ce'^S^Sy) can turn 
out to be relatively strongly negative and come into play already at e = 1. A more 
plausible explanation is that for the SV2 mode the employed asymptotic expansions 
of the operator of linearisation L, its eigenvectors and eigenvalues are valid for much 
smaller a and e than for other instabilities, because the values of Sx and Sy at the 
SV2 boundary, maximising the eigenvalue, are relatively large. Note, that the same 
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asymptotic expansion was employed in other analytic studies of the problem, cited 
in the Introduction. 

9 Conclusion 

We have presented a complete analytical study of stability of rolls near the onset of 
convection to perturbations, which are doubly periodic in horizontal directions. In 
all earlier studies only instabilities of rolls to certain classes of perturbations were 
shown. 

In pursuit of this goal, we have, first, shown that without any loss of generality 
any instability mode is responsible for either the small-angle, or finite-angle insta- 
bility. Second, for the small-angle instability modes we have derived inequalities 
determining regions of stability of rolls. The problem involves four small param- 
eters; while deriving the instability conditions we have considered all asymptotic 
relations between the small parameters. Finally, we have calculated boundaries for 
the finite-angle instability; it turns out that consideration of finite-angle instability 
modes does not modify the region of stability of rolls. 

In our analysis only the asymptotically largest terms have been taken into ac- 
count. A question often arises, whether enough terms of asymptotic expansions have 
been calculated at various intermediate stages. In 1^5^ orders of the omitted terms 
in the matrix are given, implying that the omitted terms in expressions fl66|) and 
f l68|) . used here to analyse stability, are irrelevant sufficiently close to onset. 

This small-angle instability of rolls was studied before, and the SV, OSV and ZZ 
instability boundaries found here coincide with the earlier results. Our novel results 
concerning the stability boundaries include the following ones: We have examined 
the dependence of the SV2 boundary on P for 0.543 < P < 0.782. For decreasing 
P, the boundary of the region of stable rolls on the {k, R) plane moves to the left 
away from the point {ks, -R*). We have established the asymptotics of the maximum 
growth rates and the associated eigenmodes (see (146|) ) considering exhaustively dif- 
ferent relations between a and e. We have derived asymptotic equations describing 
the regions of the instabilities and estimated remainders in these equations. 

The approach that we have followed can be applied to study instabilities of 
stripe patterns with respect to large-scale perturbations in a generic system, where 
a large-scale neutral mode exists. Existence of the invariant subspace relies only 
on the structure of equations of convection, where the linear part preserves wave 
vectors and nonlinearity is of the second order. Equations (( 16T]) and ( l63l) ) defining 
the entries of matrix A are general, they remain valid for any other system defined by 
arbitrary mappings Lj. Stability is analysed by examining the inequalities (123|) and 
( 124|) . This analysis is, perhaps, the most difficult part. It may change significantly 
for other systems with different asymptotics of the entries of the matrix A, resulting 
in different asymptotics involved in the inequalities defining instability regions. 
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A Calculation of the matrix A 

In this Appendix we calculate in the leading order the entries of the matrix A of 
the restriction of L on the invariant subspace spanned by Wj, j = 1,2,3, and 
expressions for det^ and S'(^)tr^ — det^ used to deduce the stability properties 
of rolls. The matrix, the operator and the basis are expanded in a power series in e, 
whose coefficients depend on small parameters a, 6x and 6y. Note that by virtue of 
fl2Tl) . fl22l) the action of the mapping (S^jSy) — {—S^, —Sy) amounts to permutation 
of indices W2,o ^ Ws^o- Consequently, 

Au{Sx, Sy) = y4i3(-(5^, -5y), v42i(5x', Sy) = y43i(-5^., -5y), 

(58) 

A23{Sx, Sy) = A32{-Sx, -Sy), A22{Sx, Sy) = Assi^-Sx, -Sy). 

Vector fields Wj^o (|20 i) -( l22ll . representing terms of order zero in e in the series 
(fTTI) . are eigenfunctions of Lq: 

LoW,,o = A,-oW,-o, (59) 

hence Ajjo = Xjfi and Aij^ = for i ^ j. The following relations were established 
in OP2008: 

Ai,o = -P(52 + 52), 

A^. = -C4((4 - A;2)2 + 4ak{kl - e)) + 0((4 - k^)\ a^k^ - k^)), j = 2, 3, 

C4 = 3P{4ge)-\ 9 = ^((27r2 - k'')k-^ + 3P). 

Here and below in this Appendix, plus is assumed in place of ± for j = 2, and minus 
for j = 3. 

The second and third terms of the series (|T3ll are 

Li(v, e) = (Uf™ X (V X v) + V X (V X Uf"), -(Uf" • V)^ - (v ■ V)Uf"'P), 

L2(v, 6) = (U^°" X (V X v) + V X (V X U^°") + POe,, -(U^°" ■ V)^ - (v ■ V)U2'"P). 

(60) 

The e order entries of the matrix are calculated from the relation 

3 

LoW,- 1 + Li W,- = A,, oW,, 1 + J2 ^..,1 W.,0. (61) 

i=l 
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Since the operator Lq is self-adjoint with respect to the scalar product 

(wi, W2) = wf- ■ wf " + Pi?,wf ■ w^'^'P, (62) 
the scalar product of (EI]) with Wj yields 

Aj,i = (W,,o, W,,o)-'(i.iW,,o,W,,o), 
which gives the 0(e) terms of the matrix 

Aji^i = ±\kb5y + bC2SJy + 0{S^, aS'^), 

Ay, 1 = {51 + 5l)-\-h5Jy7i\2e)-^ ± h'K\Ak^rHy{?,5l - 51) + 0{5\ a5^)), 

where j =2,3 and C2 = {Pk"^ + vi"^)(-P + l)^^(7i"^ + k"^)'^. The remaining entries 
Aij-i vanish. We use the notation 0(5^) = ©(^^^^ 5^5^ ""). 

Approximations to Wj 1 for j = 2, 3 are also found from (ICTj) . The O(e^) entries 
are calculated using the equation 

3 3 

LoW,- 2 + LiW,- 1 + LaW,- = A„ qW,- 2 + J] A.j,iW,,i + ^ 2W,,o; (63) 

i=l i=l 

the non-vanishing terms are 

Ajj^2 = -C3 ± H.iknY'S. + 0(^^ a5), 
Aji,2 = -C3 ± if2(A;vr2)-i5, + 0(5^ «5), 

where i, j = 2, 3, C3 = PP(P + l)-ia-2 and Hi + H2 = 4P(27(P + 1))-^ (in what 
follows only this sum is important). 
Finally, entries of the matrix A are 

An = -P{5l + 5'^) + 0{e'5% 

A21 = ^kbe5y + ehC25x5y + 0(e5^, ea5'^, e^), 

A31 = --khe5y + ehC25Jy + 0{e5^, ea5'^, e^), 

A12 = e(-^5.5, + ^5,(35^-5J))(5^. + 5,^)-l + 0(.5^.a5,5=^), 

A22 = -e^Cs + e^^Hi5,-C^(^{kl-ef + Aak{kl-e)'^ 

+0{{kl - P)3, a\kl - P), eV, 5^), (64) 
A32 = -e'Cs + e'^H25, + 0{e'5^e'), 

Ai3 = e[-^^^5Jy-^^^5y{35l-5l))i5l + 5'y)-' + 0{e5^ea5,e 
A23 = -e'Cs-6^^H25, + 0{e'5',e'), 
A33 = -e'C,-e^-^Hi5^-C,(^{e_-ky + Aak{e_-e)'^ 
+0{{kl - ef,a\kl - P),eH^,e^). 
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In the new basis Wi = Wi, W2 = W2 + W3, W3 = W2 — W3 the matrix of 
the operator L is 

In = -P{5l + 5l) + 0{eH'), 

Au = -ebSJy—idl + dly' + OieS^e^), 

A32 = -AC,k5^{5l + 2ak)+e'-^iH,-H2)5^ + 0{aX,Sl5J^y,e^5\e''5), 
^13 = e—6y{36l-6l){6l + 6l)-' + 0{e6\ea6,e'), 

I23 = -^C,k6,{6l + 2ak)+e^^iH, + H2)6, + 0{a%,6l6j'y,e'6^e^6), 
A33 = -CiiAk^6l + 5^ + Aak6^y) + Oiai5l + 6^),5i6le^5^,e^). 

(65) 

When calculating ( 165|) with the use of flM|) , relations (!58|) were employed to estimate 
the omitted terms. 

From (1651) we obtain 

det A = doi5l + do25i5l + c^os^^^^ + do^df 
+e\5l + 5l)-\d^^5l + ^125^5^ + ^^^^2jA + ^^^^8) 

+ 5l)-\d2,5l5l + ^22^,^) + a{d,,5X + 42^^) 

+«2(rf^^54 + ^^^5252 + ^^^^6) + ^^2(^^^^2^2(52 + ^2)-l + ^^^^4) 

+0{5\5l + 5,^)^ eH\ e'5\ a5l{5l + 5,^)^ aH\5l + 5,^)^ ae^^,^^^^ + ^4))^ 

(66) 
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where 



"01 — ~ / n , 1 \o ) "02 — "01 5 "03 — n / r-, , T NO ; "04 



(P+l)2'^"^ -ui, -u. 7r2(P+l)2' ^'^^ 7r4(P+l)2^ 

16P3 _ 32P(-3P2 + 5P + 1) 

" ~ "97r2(P + 1)2' - 277r2(P + l)2 ' 

_ 16P(-P2 + 2P + 2) _ -8P(P2 + 2P + 2) 
~ 97r2(P + l)2 ' "^1^ ~ 97r4(P + l)2 ' 

224P _ 32P 64P3A; 



(67) 



2437r4(P + l)' 8l7r4(P + l)' 7r2(P + l)2' 

32P3A; _ 16P3 J _ J J 32P3 

"32 — . / -n , 1 \9 ' '^41 — / „ , 1 \9 ' '^42 — "41, "43 



^4(p + l)2'-- (P+l)2'-^'^^ -^-^^ 7r2(P+l)2' 

128Pfc _ 32P(P2 + 2P + 2)fc 

~ 97r2(P + l)' '2 ~ 97r4(P+l)2 ' 

ti AS{A) - det ^ = eoi^^ + eo2(5^5^ + 603^^^^ + 604^^ 

+ 0(5^(52 + 54)2^ ^2^6^ ^454^ ^6^2^ ^^4(^2 + ^4)^ ^^2^2(^2 + ^4))^ 

where 

8P3(P + 5)2 _ 16P3(P + 5) _ 8P=' _ 4P3 

Cqi — / n , 1 \ ■:! — ' ^02 — 77r~Tv> — ' ^03 — ~ / -n, . 1 \ ' ^04 



(6^ 



(P+l)3 ' (P+l)2 ' (P+1)' 7r2(P+l)^ 

4p3(p2 ^ ;^gp ^ gg^) 4P(-9P^ - 54P3 - 99P2 + 20P + 2) 



en = ^ n,^ , , ei2 



97r2(P + l)3 ' 277r2(P+l)3 

4P(3P3 + 3P2 + 22P + 26) _ 4p2(p2 + 2P + 4) 
~ 277r2(P + l)2 ' ""i^ ~ 97r2(P+l) ' 

16P3(P + 5) _ 32P(P + 3)(3P2 + 2P + 2) 
~ ~8l7r4(P + l)3' ~ 2737r4(P + 1)3 ' 

16P3 _ 512P3 128PA; 

^23 — ~7r^ — /I / n , 1 \9 ' ^3 — ^77r~TT^' ^4 



8l7r4(P + l)2' " (P + l)3' " 97r2(P + l)' 

(69) 

B A necessary and sufficient condition for exis- 
tence of eigenvalues of a 3 x 3 matrix, which 
have positive real parts 

In this Appendix we show that instead of direct calculation of eigenvalues, in order 
to study stability of rolls it suffices to check, whether any of the inequalities (l23l) or 
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(|21I) is satisfied for some 5^ and 5y. 

We start by exposition of three lemmas about eigenvalues of a 3 x 3 matrix. 

Lemma 1. Let ^ be a 3 x 3 matrix with real entries. Denote its eigenvalues by A^, 
i = 1, 2, 3, and the sum of the second order minors by S{A): 

S{A) = AnA22 - AuA2i + AuA^^ - Ai^A^^ + A22A3S - 

Consider the following statements: 

SI: max(ReAj) > 

l<j<3 

S2: det^ > 

S3: S(^)tr^-det^> 

Then 

(i) S2 ^ SI and (ii) S3 SI. 

Proof, (z) follows from the identity det^ = A1A2A3 (consider separately two cases: 
all eigenvalues are real, or two of them are complex conjugate). 
Similarly, (ii) follows from the identities 

S{A) = A1A2 + A1A3 + A2A3, tr ^ = Ai + A2 + A3 

(again, consider separately the two cases). 

If in the condition of the Lemma the signs ">" are replaced by ">", the state- 
ments remain true. The modified lemma is referred to as Lemma 1'. 

Lemma 2. Let ^(x) be a 3 x 3 matrix with real entries continuously depending on 
X G r2 C R", where is a connected domain in R". Denote by Aj(x), i = 1, 2, 3 the 
eigenvalues of ^(x). Suppose 

3 xq G i7 maxReAj(xo) < 0, 

l<j<3 

det ^ 7^ Vx G Q, 
S{A)tTA-detA^O, Vx g fi. 

Then 

max ReAj(x) < Vx G 

l<i<3 

Proof. Suppose there exists xi G such that an eigenvalue of ^(x), say, Ai, has a 
positive real part. A curve in Q connects xq and Xi. The eigenvalue Ai is a continuous 
function on this curve (because roots of the cubic equation det(^ — XI) = are 
continuous functions of its coefficients.) Since ReAi has different signs at xq and 
xi, there exists a point x on the curve such that ReAi(x) = 0. If Ai(x) = 0, then 
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det^(x) = 0, and if Ai(x) is imaginary, then S{A{'k))tT A{Si) — det^(x) = 0, since 
in this case 

5(^)tr^-det^ = ReAi(2A^ + 4ReAiA3 + 2AiAi) = 2ReAi((A3 + ReAi)^ + (IniAi)^). 

Thus a contradiction with the statement of the lemma is obtained and the lemma 
is proved. 

Lemma 3. Let ^(x) and Q be the same as in the statement of Lemma 2. Assume 

3 xq G ^7 maxReAj(xo) < 0, 

l<j<3 

det ^ 7^ Vx G n, (70) 
3 xi G 1^ 5(^(xi))tr^(xi) - det^(xi) > 0. 

Then 

3 X2 G ReAi(x2) > 0, ImAi(x2) ^ 0. 

Proof. A curve in Q connects xq and Xi. Let ^ be a parameter along this curve, 
^ = at xo and ,^ = 1 at Xi. There exist ^o, < < 1, and ^i, < G < 1; such 
that 

SiAmtiAi^o) - det^(eo) = 0, (71) 
SiAiO)tTAiO - detAiO > V e G (^0,^1], 
S{A{0)tTAiO - detAiO < V e G [0,U 
Suppose all eigenvalues of ^(^0) are real. By Lemma 1' all of them are non- 
positive in [0,^o]- Due to (!70l) they do not vanish at ^0, hence they are strictly 
negative at ^o- But then (17T|) can not be satisfied. Hence the assumption that all 
eigenvalues are real is wrong. 

Let Ai and A2 = Ai be a pair of complex eigenvalues and A3 be real. By continuity, 
there exists ^2, (.2 > C,o, such that ImAi(,^) 7^ for all C, G [^05'C2]- The expression 

S{A)tTA- det A = 2ReAi((A3 + ReAi)^ + (ImAi)^) 

is positive only if ReAi is positive. Consequently, ReAi(,^3) > and ImAi(,^3) 7^ 0. 
The lemma is proved. 

The Lemmas are applied to investigate stability of rolls. 

Let A be the matrix calculated in Appendix A. Assume {6x, Sy) is the parameter 
X employed in Lemma 2, e, a and P being fixed. If (|23l) or (^M is satisfied for some 
{Sx,Sy), Lemma 1 implies existence of an eigenvalue with positive real part. 

Suppose fl23|) and are not satisfied for any 6^ and 6y. For sufficiently large 
6x (5^ ^ and 5^ a^) the matrix has three real negative eigenvalues. Let fl be 
R2 with the origin excluded. Conditions of Lemma 2 are satisfied, hence for any 
{6^, 5y) all the eigenvalues have negative real parts. 

Note that fl2^ implies that the matrix A has a real eigenvalue with a positive real 
part, while (12^ does not guarantee that there exist a pair of complex eigenvalues 
with a positive real part. However, assume in addition that det^ < for all {5xj Sy) 
(or for all {Sx, 6y) in a connected region Q, where conditions of Lemma 3 are satisfied). 
Then by Lemma 3 there exists a point {6x, Sy) where A has a complex eigenvalue 
with a positive real part. 
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C A bound for a < 0, P > Pi, ^ 

We prove here that under the conditions, stated in the title of the Appendix, 

Do + ^'^1 + £'42<(5' + 5^)"' + + a' A + 0Le^d^2^\ (72) 

is negative for all 5^ and by. 

To begin with, note that the terms involving a are asymptotically small and 
therefore are neglected. 

Assume 5^ ^ 5^. In the leading order ( 1721) is 

5j(do25^ + ^^(e^dia + rfos^J) + e'ci22 + ^'^u^J + ^045;)- (73) 
This is a quadratic polynomial in 5^, which admits a maximum at 

.9 c^os^^l + £^o?i 



^ 2(io2 

The maximum is 



(74) 



^£^(2^/03^^13 - 4rfo2t^l4)5j + e\d\^ - Mo2d22] 



4rfo2 

where both expressions in the brackets are negative for P > Pi. 

Now assume ~ Sy or ':$> Sy. In the leading order (!72|) is equal to 

(rf0l5^ + ci025X')+^'(^ll'^' + rfl25^5j + cil35^^:)(5.' + 5j)-^+^'«J^ (75) 

This quadratic polynomial in e"^ can take positive values only if 

du6l + duStSl + duSl6^ > 0. (76) 
If this is satisfied, the maximum (in e:^) of fl75]) is 
- {duS', + di25l5l + di^5l5l)\Ad225l)-' + doi^,' + (^oi + d.^Kdl + d^^^X- (77) 



In view of the inequalities (1761) . (in < and di^ > du for P > P] 



Note that (igi = c?02 and c/^g < 4(i22'io2 for P > Pi, hence fl77|) is always negative for 
P > Pi. 

D The large P limit 

In this Appendix we calculate an instability boundary, which is important for a < 
and large P. In the limit of large P the coefficient /2 in (j35ll vanishes and instability 
occurring for ~ e'^ may compete with the instability defined by (l35i) . 
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Suppose ~ e^. Then detv4 > for a < 0, if either 5^ > 5^ or 5^ > 5^ (due 
to the presence of the term ae^D^. 

Suppose 5^ ^ 51. Represent det^ (166|) (it is simpler to calculate this directly 
from (1M|) ) as 

det^ = 

- 5,^(2^2^3 + C,0{e%'^^ + (78) 

2 1,2 2 

+ (5252(_PC4A;2(£2C3 + C4O + 4P«'A;2 + ^^), (80) 

where 

For = the determinant is given by (1751) . Considering 

as a quadratic polynomial in 5^, we find that fl78l) is positive (we are interested in 
large P's, and for them the instability boundary is defined by the second term in 
(Hi ) for 

Note that 

lim /s = 00. (82) 

P— >-oo 

For large P and e not satisfying f ET]) . the contributions to det^ from (jTHD and 
(1801) are negative (the proof is omitted) . Hence if det ^ < for 6^ = and all 6y , it 
remains negative for all 6x and 6y. Thus (IHTj) is indeed a boundary for stability of 
rolls. 

For 6l ^ 5y the instability boundary is that of the Eckhaus instability 

(or in a more familiar form R — < 3{Rc{k) — R'^) ), which is below the boundary 
defined by 
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Figure 1. The area of stable rolls (shaded) on the {k,R) plane found numerically 
(see section 8) and the instability boundaries found analytically for P = 0.6 (a) 
P = 0.7 (b), P = 2 (c) and P = 7 (d), P = 20 (e) and P = 50 (f). Solid line denotes 
the onset of convection, dashed lines instability boundaries SV and OSV defined by 
fl32l) and fl35|) . dotted line the ZZ boundary fl36|) and dashed-dotted line the SV2 
boundary fHTl) . fH2]) . Stars mark the points where the SV2 boundary is obtained by 
interpolation. Horizontal axis: k, vertical axis: R. 
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Figure 2. The dominant growth rate (sohd hne, left vertical axis) and the values of 
Sx and 5y (dotted and dashed lines, respectively, right axis), where the maximum is 
achieved, versus the Rayleigh number (horizontal axis) for P = 2 and A; = 2.15 (the 
respective crossection is shown by a thin vertical line on fig. Ic). 
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